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We study the Josephson effect between atomic Bose- 
Einstein condensates. By drawing on an electrostatic anal- 
ogy, we derive a semiclassical functional expression for the 
three-dimensional Josephson coupling energy in terms of the 
condensate density. Estimates of the capacitive energy and of 
the Josephson plasma frequency are also given. The effect of 
dissipation due to the incoherent exchange of normal atoms 
is analysed. We conclude that coherent Josephson dynamics 
may already be observable in current experimental systems. 

Recently jfj-^ , it has become possible to cool a macro- 
scopic number 10^ — 10^) of magnetically confined, 
spin-polarized atoms down to temperatures of order 100 
nK while maintaining densities sufficiently high (10^^ — 
lO^^cm^"^) to permit the onset of Bose-Einstein con- 
densation (BEG). From the ensuing theoretical work, it 
has been concluded that these Bose condensed atomic 
gases behave very differently from the ideal noninteract- 
ing gases, which yields the prospect of potentially dis- 
playing a rich phenomenology that might include vortex 
states and the Josephson effect . 

In this article we present a theoretical study of the 
Josephson dynamics between two atomic baths that have 
undergone BEG. The Josephson effect results from a col- 
lective mode of two weakly connected systems between 
which a macroscopic fraction of particles can tunnel with 
identical probability amplitude. Going beyond previ- 
ously proposed one-dimensional and dissipation free 
10,^ models, we present here a three-dimensional study of 
the Josephson effect between Bose condensates and esti- 
mate the effect of damping. We calculate the Josephson 
coupling energy, the capacitive energy (which accounts 
for quantum fluctuations of the phase) , and the frequency 
of the Josephson plasma oscillation. Our main conclusion 
is that still lower temperatures than those achieved up to 
date are needed for a clear realization of the Josephson 
effect. 

The collective dynamics of a Bose condensate at zero 
temperature is described by its macroscopic wave func- 
tion ^'(r, i). If this is factorized as 5" = .ypexp(i</?), the 
standard energy functional can be written as 



of the phase across the condensate near equilibrium in a 
single spherical harmonic well [see Eq. (|l|)]. 
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{g — 2TTh^a/m), and the corresponding Hamilton equa- 
tions lead to the Gross-Pitaveskii equations |||,^. Ne- 
glecting depletion [|l^, the normalization can be taken 
as / drp = N, being N the total number of atoms. 

At sufficiently low temperatures the phase within one 
well can be regarded as uniform. This can be easily seen 
if one estimates the energy of a one-radian fluctuation 
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where R = ao(15aAf/ao)^/^ is the cloud radius esti- 
mated within the Thomas- Fermi approximation jj], ujq 
is the harmonic oscillator frequency within the well, and 
Go — {h/mujoY^'^ is the oscillator length. For the last 
approximate equality we have used typical parameters 
a = 5 nm and ao = lO"** cm. The characteristic temper- 
ature of such an oscillation can be as big as 10 /iK for 
the MIT (1996) experiment §. The estimate (|) sug- 
gests that, in the interestitial region between two wells, 
where p(r) decreases appreciably, spatial phase variations 
are less costly and thus easier to activate at low temper- 
atures. Here lies the essence of the low energy Josephson 
dynamics which we study below. 

Let us consider two weakly connected condensates 1 
and 2. We assume that the condensates are confined 
within spherical harmonic wells of the same frequency 
ujQ. First we wish to analyse the semiclassical dynamics. 
It is, of course, well known that the relative number SN 
and the relative phase x of condensate in the two 
wells may be treated as canonically conjugate variables. 
Nevertheless, by introducing a sufficiently coarse-grained 
average n of the number in (say) well 1 we may treat n 
and X simultaneously well-defined, and write for the 
wave function the ansatz 

*(r, t) cx *i(r; n{t)) + e*>^(*)^'2(r; N - n{t)), (3) 

where ^'^(r; n) is the (real) equilibrium wave function for 
the isolated well i containing n bosons. It is straightfor- 
ward to show that, to lowest order in the overlap integrals 
/ 5'i'I'2, the energy fmictional for 5* takes the form 

H{SN,x) ^ Eb{SN) + Ej{SN) (1 - cosx) , (4) 

where Eb (SN) is the bulk energy of the two isolated wells 
with 5N = N/2~n transferred atoms, and Ej{5N) is the 
Josephson coupling energy. Eb (SN) may be expanded as 



Eb{SN)^Eb{0) 
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Within the Thomas- Fermi approximation, ~ N'^/^ [|j, 
so that the ratio between the third and second terms in 
the expansion is 0.32 SN'^/N'^, which means that the last 
term can be neglected in a wide range of situations. To 
avoid complications stemming from possible resonances 
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between Josephson oscillations (see below) and intrawell 
excitations, we require yu(iV/2 + SN) - ijl{N/2 - 6N) < 
hujQ, where we use the result that the first normal mode 
of a spherical well lies approximately at Hloq above the 
ground state [^,0. This condition is realized when 
SN/N ^ 4.6A^ for typical parameters. This may 
seem an important restriction on SN. However, the frac- 
tion of transferred atoms can be as big as 10% for the 
JILA experiment , but it has to be < 2% for the MIT 
experiment jsj. 

The expression for the coupling energy Ej{dN) in 
terms of and ^'2 which can be derived from Eqs. (||- 
^ is rather complicated and difficult to handle. A much 
simpler expression can be obtained if one notes that the 
coupling energy must come entirely from the term in 
Eq.([l]), which in turn can be approximated as 



/ dr|lp(r,t)|V^(r,i)pEEi?M, 
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where the integration extends over the the region exterior 
to the condensates (the results are quite independent on 
the precise location of the condensate borders). As ar- 
gued before, the phase within the condensates can be 
assumed to be practically uniform. In the limit of very 
small phase difference, x ^ 1; the phase term in Eq. (Q) 
can be approximated as Ejx^ /2. The only way for Eq. 
(^ to have such a dependence on the total phase differ- 
ence X Slid not on the details of V(/? is that the condi- 
tion SE[Lp]/S(p{r,t) = is satisfied. One may note that, 
except for trivial factors, this is the electrostatic equa- 
tion for the electric displacement vector D = pV(p in 
a medium with a nonuniform dielectric constant p{r, t) . 
Boundary conditions for ip{v, t) are given by its value 
at the borders of each condensate, which act as con- 
ductors in this analogy. We have a system of two con- 
ductors held at a potential difference x and a dielectric 
medium surrounding them. Then E[(p\ is essentially the 
energy of this capacitor. Potential theory tells us that 
E[ip\ ~ (fi^ /2m)C[p\x^, which imphes 



Ej = -C[p], 
m 
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where C[p] is the mutual capacitance of the two conduc- 
tors in the presence of such a dielectric jl^ . 

We take the origin of coordinates in the middle point 
of the double well configuration, taking the z direction 
along the line which joins the two minima. Standard 
variational arguments can be invoked to prove that, for 
parallel plate boundary conditions. 
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The lower bound is obtained by removing the positive 
term (dip/dx)"^ + {dip/dy)'^ from the energy functional 



(^, while the upper bound is derived by taking (p in- 
dependent of x, y. If Vext depends weakly on x, y in 
the region that controls the capacitance, we can write 
p{x,y,z) ~ p{0,0,z) and then the two bounds become 
approximately equal to 
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where A is an effective area. 

In order to proceed further, we need an estimate of p = 
I^Pp in the region of interest. For two identical wells in 
equilibrium, the ground state wave function is symmetric 
in z. Neglecting the dependence on x, y, and within the 
WKB approximation , we can write 
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where p{z) — [2m(14xt(0, 0, z) — p)]^^^ and B is a con- 
stant to be determined later. Introducing (^o|) into (H) 
we obtain 
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where 5*0 = J^p{z)dz/h. For standard (quartic) barriers, 
Sq ~ 27ra(Vb — p)/huJo, where Vq is the barrier height 
along the line x — y ^ 0, and a is of order unity {a = 
8/37r for < Fq and a = l/%/2 for Vq - p <. Vq). 

The coefficient B must be calculated by properly 
matching ( p^ with good approximate solutions near 
the border of each well. If 6*0 ^ 1, then Eq. ( p^ ) 
can be matched with the solution given in Ref. i7i for 
the region R:s>z~R:s>d, d = uq {ao/2R)^'^ be- 
ing the distance from the classical radius R where the 
Thomas-Fermi approximation begins to fail. The result 
is B ~ ue-^''/2(n/87rad3)i/2^ with u ~ 0.397 0. To es- 
timate the effective area A we note |^] that the form of 
the order parameter has a universal dependence on the 
position near the boundary of the Thomas-Fermi zone, 
namely ^'(r) ~ [(r — R)/d\ /d^/Sna, where </> is a uni- 
versal function (in that it does not depend on the confin- 
ing potential) which varies appreciably within a scale of 
unity. Therefore p{x, y, z) varies transversally on a scale 
such that {\J x"^ + + i?^ — R)/d ^ 1. These considera- 
tions yield an estimate oi A = 2^/'^w(ao/i?)^/'^7ri?^, where 
-y ~ 1. The result is 
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We note that, because the prefactor in has the same 
N^/^ explicit dependence as the critical temperature 
for a system of free confined bosons (which is approx- 
imately equal to the critical temperature of interacting 
bosons |12[| ) , one could expect a simple relation between 
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the two magnitudes. Knowing that fc^Tc ~ N^^^Hluq 
, then (jlJ) can be rewritten as (for typical cases) 
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Also within the Thomas-Fermi approximation, a sim- 
ple expression can be obtained for Eb{SN) in Eq.(^, 
which, up to quadratic order in SN, can be written 
EeiSN) - EeiO) = EcSN^/2, where Ec = 2/i'(iV/2) 
is the capacitive energy. From the result of Ref. Q] for 
/i(iV/2), we obtain 
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Collecting terms, the Hamiltonian (0) can be written 
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£;j(l -cosx), 
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giving the well-known equivalence to a pendulum [ [L3[ . 
The frequency of small oscillations (the Josephson plas- 
mon), ujjp — y'EjEc/h, turns out to be 
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which is typically a fraction (not necessarily very small) 
of the confining potential. As usually a;o/27r ~ 10 — 100 
Hz, we conclude that, ujjp < 10 Hz. 

The ratio Ej/Ec is a good measure of the classical 
character of the relative phase x- From ( p^ ) and (14), 
we find 
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By varying 5*0, the system can be driven from the clas- 
sical regime (Ej ^ Ec) to the strong quantum limit 
(Ej <C Ec)- However, in the latter case, quantum fluc- 
tuations are only important if we operate at ultralow 
temperatures ksT < Ec- 

The Hamiltonian (|l^) describes the dynamics of a con- 
servative system. In real life, however, we should expect 
a certain amount of damping. The most obvious source 
of such damping is the incoherent exchange of normal 
atoms, and a quantitative discussion requires a general- 
ization of our results to nonzero temperature. Because 
the spatial scale of the normal component is quite dif- 
ferent, in a harmonic trap, from that of the condensate, 
this generalization is much less trivial than for the case 
of a junction linking two homogeneous superconductors 
or superfluids, and we shall not attempt a quantitative 
discussion here. However, we shall give two qualitative 
arguments, based on two very different physical assump- 
tions (corresponding essentially to the high and low bar- 
rier limits), to the effect that the damping associated 
with the normal component will overdamp the Joseph- 
son behavior and hence make it in practice unobservable 



down to temperatures of the order of huJo/kB- In both 
cases, we find that normal atoms give an Ohmic contri- 
bution to the current. 



In = -GSfi, 
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where 6 fi is the chemical potential difference (for simplic- 
ity, we assume equilibrium within each well). As a result, 
the first Josephson equation is modified to read 



—dN^ — sinx + /„. 
at n 
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In the high barrier limit (Vb ^ ksT), the basic order- 
of-magnitude assumption is that any power-law factors 
occurring are negligible compared to the relevant WKB 
or Arrhenius-Kramers exponentials. At temperatures 
T < Tc, the normal component in each well will be ap- 
preciable and will be distributed over an energy range 
- keT > huJo, since @ kBTc/tuvo ~ N'^^^ > 1. We 
will assume that fc^T is, nevertheless, small compared to 
Vi = Vq — ^0, with /ig = /i(iV/2). The situation we have 
is Vq ^ ksT ^ fiQ ^ ^^0- Under these conditions the 
typical spacing of the one-particle energy levels at ener- 
gies ksT is small cornpared to fiojQ, and we can ignore 
"level-crossing" effects ||lj] and treat the tunneling of un- 
condensed particles as incoherent. At high T, the total 
rate of crossing is then given by a standard Arrhenius- 
Kramers formula, P„ = (KijJo/27r) exp(— T/i/ZcsT), where 
usually |l^ K ^ 1. We have, approximately, G ~ 
PnNn/ksT where iV„ is the number of normal particles. 
The degree of damping can be estimated by comparing 
( p8| ) with the Josephson supercurrent in a small oscilla- 
tion of amplitude (5/x, namely. 
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From dlj), we estimate ojjp ' e-^«i^ujo {Nq - N 

is the number of condensed particles), valid for typi- 
cal parameters at T = 0. Finite T corrections should 
not change the qualitative conclusions. We find Is /In ^ 
2TrNl^^^{nLOo/kBT)^exp{Vi/kBT-So/2), since N^/N - 
{T/Tc)^- We conclude that, in the high-barrier limit. 
Is I In ^ 1; hence the motion of the equivalent pendulum 
is overdamped, at least down to temperatures of the or- 
der of hivo/ks- For temperatures below this our estimate 
fails for several reasons, not least because the density of 
normal component falls exponentially rather than as T^. 

The estimate presented above is, however, quite irrel- 
evant for today's experiments, because the high barrier 
condition requires Sq/tt ^ 2fio/fkuo ~ 6 — 100 (we take 
a ~ 1) for typical situations. On the other hand, the 
prefactor in Ej/h can be estimated as 10* — 10^ Hz, and 
thus the critical current itself would be so small as to be 
completely unobservable. A more relevant regime is that 
in which the chemical potential lies near the top of the 
barrier, Vi fiujQ/2, so that S'o ~ 1 — 5. In this regime, 
we can still expect the WKB formula for Ej derived 
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above to yield a reasonable approximation. However, if 
ksT huJo/2 ~ Vi, the thermal cloud lies mostly above 
the barrier and a radically different approach is needed to 
study its transport properties. For simplicity, we intro- 
duce the drastic approximation that particles impinging 
on the barrier with energy E are transmitted with proba- 
bility one if £■ > Vb, and zero ii E < Vq. Then, the flow of 
normal atoms due to a fluctuation in Sfj, is only limited by 
the "contact resistance", a concept taken from ballistic 
transport in nanostructures [p^ . Adapting standard ar- 
guments to the case of bosons, we may write G = Nch/h, 
where A^ch is an effective number of available transmis- 
sive channels. Within a continuum approximation, and 
assuming that only transverse channels with a minimum 
energy between fiQ and /io + ksT are populated, we find 
(taking Sfj, <C ksT) Nch = mAnkBT /2'kIi' , where An 
is a mean transverse contact area seen by normal parti- 
cles ll^. Approximating An ^ 2'KkBT /mui'^, we find 
A'ch ~ {ksT /Tlloq)^ , and hence 



Is/In - 2^Ny^\hujQ/kBTfe-^"'^, 
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which, interestingly, is formally equivalent to the high 
barrier result (since there 5*0/2 ~ ttoViITiloq ^ Vi/kBT). 
Taking kuT = lO^iwo, we find h/In ^ 0.38 for A^ - 10*, 
and 3.3 for A^ ~ 10^, if ~ 5. The equivalent numbers 
for 5*0 ~ 1 are 2.8 and 24. We tentatively conclude that 
effects of coherent Josephson dynamics can be observed 
in today's atomic Bose condensates if the barrier is low. 
Underdamped dynamics can be further favored by de- 
creasing T and increasing A^. It is important to note, 
however, that some aspects of the Josephson behavior 
can be observed in the overdamped regime, provided that 



thermal fiuctuations are unimportant [ksT/Ej 



< 



1), 



[ |18| a not very stringent condition in the low barrier limit 
[seeEq. (H)]. 

In summary, a systematic study of the Josephson ef- 
fect between two weakly connected atomic Bose-Einstein 
baths has been presented. We have derived a novel, 
three-dimensional functional expression for the Joseph- 
son coupling energy in terms of the condensate density. 
The capacitive energy and the frequency of the Josephson 
plasma oscillation have also been calculated within the 
WKB approximation. The effect of damping due to the 
incoherent exchange of normal atoms has been estimated 
in the limits of high and low potential barrier. Within 
the low barrier regime, we find [see Eq. (|2l|)] that weakly 
damped Josephson dynamics may already be observed in 
current experimental setups, and that coherence between 
atomic Bose condensates can be further enhanced by low- 
ering the temperature and the potential barrier, as well 
as by increasing the number of condensate particles. 
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